Draft version February 2, 2008 

Preprint typeset using LATpjX style cmulateapj v. 11/12/01 



NGC 4261 AND NGC 4697: REJUVENATED ELLIPTICAL GALAXIES 
A. Zezas, L. Hernquist, G. Fabbiano, J. Miller 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 
azezas@cfa. harvard. edu;lhernquist(Qcfa. harvard.edu; pepiOcfa. harvard. edu;jmmiller@cfa. harvard.edu 

Draft version February 2, 2008 

ABSTRACT 

Chandra images present evidence for a non-uniform spatial distribution of discrete X-ray sources in the 
elliptical galaxy NGC 4261. This non-uniform distribution is inconsistent with the optical morphology 
of NGC 4261 at greater than a 99.9% confidence level. Similar evidence is seen in one more elliptical 
galaxy (NGC 4697; 98% confidence level) out of five cases we investigated. NGC 4261 and NGC 4697 
have old stellar populations (9-15 Gyrs) and fine structure parameters of 1 and respectively, suggesting 
no recent merging activity. On the basis of simulations of galaxy interactions, we propose that the X- 
ray sources responsible for the non- uniform distribution are associated with young stellar populations, 
related to the rejuvenating fall-back of material in tidal tails onto a relaxed merger remnant, or shock 
induced star-formation along the tidal tails. 

Subject headings: galaxies: elliptical and lenticular — galaxies: individual(NGC 4261; 3C270) - 
galaxies:interactions — X-rays: galaxies — X-ray: binaries 



1. INTRODUCTION 

Numerical simulations of galaxy interactions suggest 
that mergers of spiral galaxies can lead to the formation 
of elliptical galaxies (e.g. Toomre & Toomre 1972; Barnes 
1988, 1992; Hernquist 1992, 1993). These expectations 
are supported by optical imaging observations which show 
that several elliptical galaxies exhibit shells, ripples, arcs, 
counter-rotating cores, or faint tidal tails; structures which 
are interpreted as evidence for galaxy interactions or merg- 
ers (e.g. Schweizer et al. 1990). 

The availability of high spatial resolution X-ray data al- 
lows us to investigate these indications using X-ray binary 
populations as probes of the star-formation histories of 
galaxies: High-Mass X-ray binaries (HMXBs) form much 
more efficiently than Low Mass X-ray binaries (LMXBs) 
(Kalogera & Webbink, 1998; Portegies Zwart & Verbunt, 
1996), leading to a higher number of X-ray sources per star 
in young stellar populations. Thus, using number counts 
of discrete X-ray sources, we can identify regions of recent 
or enhanced star-formation within a galaxy. For example, 
the spatial distribution of LMXBs (which trace the old 
stellar populations) in early-type galaxies is expected to 
be smooth, generally following the distribution of optical 
star-light. On the other hand, if there are any sites of re- 
cent star-formation (e.g. triggered by galaxy interactions) 
they are expected to host HMXBs, the larger numbers of 
which (compared to the LMXBs) may result in an overall 
non-uniform X-ray source spatial distribution. 

In this paper we present evidence for a non-uniform spa- 
tial distribution of X-ray sources in NGC 4697, an ellipti- 
cal galaxy without any indication for merging activity from 
optical data, and we compare these results with similar ob- 
servations of other early- type galaxies, finding a second ex- 
ample in NGC 4261. A more detailed investigation of the 
X-ray source populations, their degree of non-uniformity 
and links to other merger activity diagnostics will be pre- 
sented in a forthcoming paper. 



2. NGC 4261 

NGC 4261 is a nearby (35.1 Mpc for H G =75 km/s) 
E2 type galaxy, belonging to a poor group of galaxies 
(Garcia, 1993) but without any evidence for interaction 
with other members of the group. Optical observations 
show a 20"(^ 3.4 pc) long dust lane along the North- 
South axis of the galaxy (Martel et al. 2000 and references 
therein) and a spectacular nuclear dust disk (~ 2"; Jaffc 
et al. 1996, Ferrarese et al. 1996). Its optical isophotes 
in larger scales have a very strong boxy morphology (e.g. 
Peletier et al. 1990). Optical spectroscopy shows that 
its dominant stellar population has an age of ~ 15 Gyr 
(Trager et al. 2000). This galaxy also features two promi- 
nent radio jets (e.g. Birkinshaw & Davies 1985, Jones et 
al. 2000) emanating from an active nucleus. The mass 
of the supermassive nuclear black-hole is estimated to be 
(4.9 ± 1.0) x 1O 8 M (Ferrarese et al. 1996). 

2.1. Observations and Data Analysis 

NGC4261 was observed for 35 ks (OBSID 834; PI M. 
Birkinshaw) with the ACIS-S3 CCD on-board the Chandra 
X-ray Observatory. The spatial resolution of Chandra is 
0.5" (van Speybroeck et al. 1997) corresponding to a phys- 
ical resolution of ~ 80 pc at the distance on the galaxy. 
Although the observation was performed in 1/2 subarray 
mode in order to mitigate pile-up on the nuclear source, 
the active field of view (4.1' x 8.3') was large enough to 
cover the whole galaxy (D 2 5 = 4.07'; de Vaucouleurs et al., 
1991). The initial processing of the data and the results 
from the analysis of the nuclear spectrum are presented in 
Zezas et al. (2003). In this paper we concentrate on the 
spatial distribution of the discrete sources. 

We used the wavdetect algorithm, within the CI AO v2.2 
data analysis suite in order to search for discrete sources 
in four different energy bands: full band (0.3-7.0 keV), soft 
band (0.3-2.0 keV) and hard band (2.0-7.0 keV). The limit 
of 7.0 keV was set in order to minimize contamination by 
the particle background, while the boundary of 2.0 keV 
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was chosen in order to limit the contribution of the dif- 
fuse emission (which is mainly thermal with kT = 0.6 keV) 
only in the soft band. The source detection was performed 
in scales of 1, 2, 8 and 16 pixels and a limiting probability 
of ~ 1 chance detection over this field. We found a total 
of 62 sources within the S3 CCD. Of those, 45 sources are 
within the D25 radius of NGC 4261, 40 of which have sig- 
nificance higher than 3cx above the local background. The 
absorption corrected luminosities of the latter are above 
the Eddington limit for an 1.4 MQncutron star, ranging 
from 2.5 x 10 38 erg s _1 up to 4.4 x 10 39 erg s _1 (with two 
sources above 10 39 erg s" 1 ) assuming a power-law model 
with r = 1.7 and Galactic line-of-sight Hi column density 
(5.8 x 10 20 cm" 2 ; Stark et al. 1992). 

2.2. Spatial distribution of X-ray sources 

Figure la shows an adaptively smoothed image of 
NGC 4261 in the full band, together with the D25 ellipse 
of the galaxy and the active field of the S3 ACIS-CCD. For 
reference in Figure lb we present a DSS2 R-band image 
of the galaxy in the same scale, with the discrete sources 
marked by the best fit 3a ellipses to the spatial distribu- 
tion of their photons 1 . From the X-ray image is clear that 
the spatial distribution of the discrete sources does not fol- 
low the distribution of optical star-light, but shows instead 
a distinct asymmetric pattern. In the following discussion 
we focus on the sources within the D25 ellipse. 

In order to assess the probability that the non-uniform 
distribution may be observed by chance, we performed a 
Kolmogorov-Smirnov (KS) test to compare the distribu- 
tion of the position angles (PA) of the sources (measured 
clockwise from the North-South direction), with (a) a uni- 
form distribution and (b) the azimuthal profile of the opti- 
cal light of the galaxy. For the latter we created a sample 
of "optical sources" based on the star-light distribution 
in the DSS2 R-band image, by randomly drawing posi- 
tion angles in 2° sectors from a flat parent distribution 
(equivalent to assuming constant surface brightness within 
each 2° sector). The number of draws in each sector was 
proportional to the fraction of the total galactic light in- 
cluded in it, after removing any strong point-like optical 
sources (most probably associated with foreground stars 
or background galaxies/ AGNs). We found that both com- 
parison distributions are inconsistent with the observed 
X-ray source azimuthal distribution above the 99.9% con- 
fidence level (Fig. 2a). This result holds even for model 
distributions with as few as 10 sources. 

A search for irregularities in the azimuthal distribution 
of the PA of the sources, using the Bayesian blocks method 
(Scargle 1998; Scargle et al. 2003, in prep) also shows 
that there is an excess of sources between PA=140° and 
PA=190° at a confidence level greater than 99% (esti- 
mated from the prior of 7 = 4.5 used in the decomposition; 
however see Scargle et al. 2003, in prep 2 , for a caveat on 
this interpretation of the value of the prior). Therefore, we 
conclude that the spatial distribution of the X-ray sources 
is not the same as that of the optical light tracing the 
stellar population. 

3. OTHER EXAMPLES: NGC 4697 

1 http://asc.harvard.edu/ciao/download/doc/dctcct_html_manual/ 

2 see also http://space.mit.edu/CXC/analysis/SITAR/functions.html 



Indication for a different distribution of X-ray sources 
and optical star-light has been also seen in NGC 720 (Jel- 
tema et al. 2003), which has relatively young stellar popu- 
lations (~ 5 Gyrs; Trager et al. 2000). The X-ray sources 
in NGC 720 appear to trace out spiral arms, but perform- 
ing the same analysis as above we confirm that this result 
is not statistically significant (Jeltema et al. 2003). Based 
on the KS test we find that the X-ray source distribution is 
different from the optical distribution or a constant at the 
50% and 70% confidence levels respectively. Similarly, the 
Bayesian blocks method showed that any local enhance- 
ment of sources is significant only at the 63% confidence 
level. 

We searched the Chandra archive for deep observations 
of nearby early-type galaxies, which also did not show any 
evidence for fine structure (Schweizer & Seitzer 1992) or 
recent merging activity, in order to construct a comparison 
sample. We obtained data for NGC3379, NGC4636 and 
NGC4697, the general properties of which are presented 
in Table 1. We analyzed the data for these three galax- 
ies in the same way as for NGC 4261 (§2.1). A KS test 
between the X-ray source distribution and the compari- 
son distributions (calculated as described in §2.2) showed 
that only in the case of NGC 4697 we can rule out at 
the 98% confidence level (Fig. 2b) that they belong to 
the same population (~ 95% probability that there is a 
local enhancement of sources based on Bayesian blocks). 
In the other two cases (NGC 3379 and NGC 4636) these 
confidence levels are ~ 30% and ~ 25% respectively, not 
allowing us to draw any definitive conclusions. 



4. DISCUSSION AND CONCLUSIONS 

Our results show evidence for a non-uniform spatial dis- 
tribution of the X-ray source population in two nearby 
apparently normal elliptical galaxies. These galaxies have 
very low fine-structure parameters (S = 1 for NGC 4261 
and S = NGC4697) indicating that if they are merger 
products the merger event took place at least a few Gyrs 
ago. This is because merger simulations indicate that 
ripples, shells or strong tidal tails are not easily observ- 
able for much longer than ~ 10 9 yr after the relaxation of 
the merged system (e.g. Quinn 1984; Hcrnquist & Quinn 
1988, 1989). Moreover, faint traces of tidal tails surviving 
for much longer might not be easily detectable in optical 
observations because of their very low surface brightness 
(Mihos 1995). In fact, optical observations indicate stellar 
populations as old as 9 and 15 Gyrs in the nuclear re- 
gions of NGC 4697 and NGC 4261 respectively (Trager et 
al. 2000). 

We propose that one or more localized star-formation 
events which occurred at most a few hundred Myr ago, 
are responsible for the discrepancy between the optical 
and X-ray morphology of these two galaxies. Accord- 
ing to X-ray binary formation models, LMXBs are more 
susceptible to effects such as supernova kicks and com- 
mon envelope phases (Kalogera & Webbink, 1998) than 
HMXBs, leading to HMXB formation efficiencies 10 to 
100 times higher than for LMXBs (Portegies-Zwart & Ver- 
bunt 1996). In fig. 3 we plot the Lx/Lb ratios of star- 



Rejuvenated elliptical galaxies 



3 



forming galaxies from the Einstein survey of Shapley, Fab- 
biano & Eskridge (2002) together with the mean Lx/Lb 
ratio for the discrete X-ray sources in elliptical galaxies 
(log(L x /L B ) = -3.05 ± 0.4 3 ; Athey 2003). From this plot 
is clear that even after accounting for a 50% contribu- 
tion from the diffuse emission in star-forming galaxies (e.g. 
Zezas et al. 2001), the latter have systematically higher 
Lx/Lb ratios than the X-ray binary component of early 
type galaxies. This together with the fact that individual 
HMXBs have significantly lower X-ray to B-band flux ra- 
tios than LMXBs (van Paradijs & McClintock 1995) sup- 
ports the picture that HMXBs are forming much more 
efficiently than LMXBs. 

Recently Barnes (2003) suggested that shock-induced 
star-formation can explain the star-forming activity ob- 
served along the tidal tails of "The Mice" (Stockton 1974; 
de Grijs et al. 2003). According to this model shock waves 
developing along the tidal tails can compress the neutral 
gas and trigger star-formation. This picture is consistent 
with the tail-like distribution of the observed X-ray sources 
in NGC 4261, if the tidal tails are projected against the 
body of the merger remnant. If the age of this young stel- 
lar population is less than a few hundred Myr, it forms 
X-ray binaries much more efficiently than the populations 
in the relaxed merger remnant. On the other hand its op- 
tical emission is diluted by the optical light of the much 
stronger old population making its detection in the op- 
tical band very difficult. Depending on the strength of 
the star-formation, this may result in an overall projected 
spatial distribution of X-ray sources which is inconsistent 
with that of the star-light. This star-formation event can 
be well-approximated by an instantaneous burst, in which 
case the young X-ray source populations are expected to 
decay in a few hundred Myr after the passage of the shock. 
Therefore, we estimate that the shock should propagate 
with a speed of ~ 100 km s in order to cover the length 
of the region we observe X-ray sources (20.5 kpc), within 
the timescale of HMXB formation (~ 200 Myr). This 
speed is realistic for shock waves developing in the inter- 
stellar medium. 

Alternatively, N-body simulations show that structures 
resembling dwarf galaxies can form in the tidal tails as a 
merger completes (Barnes & Hernquist 1992, 1996). These 
structures remain bound to the body of the remnant, or- 
biting it. Because the objects formed in this manner have 
a range of binding energies in the potential well of the rem- 
nant, the characteristic time for them to fall back onto the 



remnant can be much longer than the time for the body 
of the remnant to relax (Hernquist & Spergel 1992; Mihos 
& Hernquist 1996). Thus, depending on the orbital dis- 
tribution, the remnant will eventually lose evidence for 
a merger, while tidal dwarfs continue to orbit at large 
radii. The interaction between the body of the merger 
and the condensations in the tidal tails can trigger small 
scale star-formation events and locally enrich the galaxy 
with a young stellar population. 

As was mentioned earlier, a young stellar population can 
form X-ray binaries very efficiently, while it is very difficult 
to detect it against the much brighter stellar populations of 
the merger remnant. Therefore, depending on the strength 
of the star-formation, this may result in an overall spatial 
distribution of X-ray sources which is inconsistent with 
that of the star-light. Since fall-back is estimated to occur 
for several Gyrs after relaxation, we expect that even el- 
liptical galaxies with stellar populations of ~ 10 Gyrs may 
exhibit non-uniform X-ray source distribution, if they are 
the end-points of galaxy mergers. 

Although the spatial distribution of the sources in 
NGC 4261 indicates that they are associated with the tidal 
tails, probably both mechanisms can produce populations 
of numerous, young X-ray binaries in merger remnants. 
Given the timescales of the orbits of dwarf galaxies (up to 
1 Gyr) or the velocities of the shocks along the tidal tails, 
and the lifetime of HMXBs (up to ~ 100 Myr), we would 
expect long periods during which these galaxies have uni- 
form X-ray source distributions. These periods are ex- 
pected to be longer in the more evolved systems since 
the dwarf companions with short orbits are accreted first, 
which is consistent with our finding that this phenomenon 
appears in only a few of the most evolved galaxies exam- 
ined. Studies of a more complete sample of elliptical galax- 
ies and merging galaxies in their latest stages of merging 
will allow us to further investigate this hypothesis. 
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Table 1 



Properties of the galaxies 



Galaxy 


Hubble Type 


Distance^ (Mpc) 


E* 


Notes 


NGC 4261 


E2 


32.5 


1.0 


dust, quadrupole structure (a) 










age ~ 15 Gyrs (b) 


NGC 3379 


El 


8.1 


0.0 




NGC 720 


E5 


28.0 




in group, age 5 Gyrs (b) 


NGC 4636 


EG 


15.0 




complex X-ray morphology (c) 


NGC 4697 


E6 


15.9 


0.0 


age ~ 9 Gyrs (b) 



* Galaxy distance assuming H — 75 km/s/Mpc. 

^Fine structure parameter from Schweizer &: Seitzer, 1992. 

References: (a) Colbert, Mulchacy & Zabludoff, 2001; (b) Trager et al. 2000; (c) Jones 
et al. 2002. 



Fig. 1. — Left: (a) A full band (0.3-7.0 keV) adaptively smoothed Chandra image of NGC 4261; the D25 ellipse is shown by the white 
ellipse. The white box indicates the active field of the ACIS-S3 chip. The green ellipses indicate the position of the sources in the central, 
saturated, region of the image. Right: (b) The R-band DSS2 image of the galaxy with the X-ray sources marked by the red ellipses. The 
D25 ellipse is again shown by the white ellipse. 
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Fig. 2. — (a) Comparison between the cumulative distribution of the position angles of the X-ray sources within the D25 ellipse of NGC 4261 
with the azimuthal distribution of its optical light (dashed red line) and a constant (dotted green line), (b) The same plot but for NGC 4697. 
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Fig. 3. — Histogram of the X-ray (0.35-10.- keV) to total B-band luminosity ratio of spiral galaxies later than Sc type (T> 6) from the 
Einstein sample (Shapley et al. 2001). The shaded regions show only the X-ray detections whereas the unshaded regions show both the 
detections and non-detections (upper limits). The horizontal line gives the range of Lx/Lb for the discrete source populations of elliptical 
galaxies (Athey 2003). 
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